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I. INTRODUCTION 



The QCD theory allows the existence of bound states of gluons, called glueballs, but no 
firm experimental discovery of such states has been obtained yet. An important difficulty 
is that glueball states might possibly mix strongly with nearby meson states. Nevertheless, 
the computation of pure gluon glueballs remains an interesting task. This could guide 
experimental searches and provide some calibration for more realistic models of glueballs. 

Lattice calculations are undoubtedly a powerful tool to investigate the structure of glue- 
balls. A previous study []] predicts the existence of a lot of resonances between 2 and 4 GeV. 
A recent update of this work confirms the results already obtained. 



The potential model, which is so successfu 



a possible approach to study glueballs P, [4], |5|, |6( . In a recent paper |7(, a semirelativistic 
Hamiltonian is used to compute two-gluon glueballs with masses in good agreement with 



to describe bound states of quarks, is also 



those obtained by the lattice calculations of Re 
Reta^eHesoothea^Md,—, 



This Hamiltonian, the model III in 
2| and on a one-gluon exchange (OGE) 



interaction proposed in lief. Q]. It depends only on three parameters: the strong coupling 
constant as, the string tension a, and a gluon size 7 which removes singularities in the 
short-range part of the potential. The constituent gluon mass is dynamically generated and 
it is assumed that the Casimir scaling determines the color structure of the confinement. 
These two ingredients are actually necessary to obtain a good agreement between the results 
from a potential model and from lattice calculations. 

The purpose of this paper is to check if the potential model built for two-gluon systems 
in Ref . 3| can be generalized to three-gluon systems. Compared to previous models |f| , 
our approach is characterized by some improved features: semirelativistic kinematics, more 
realistic confinement, dynamical definition of the gluon mass, coherent treatment of the 
gluon size. These points will be detailed below. The masses of the lowest negative parity 
L = glueballs are computed with a great accuracy and compared with lattice calculations 
In Sec. |TIJ the three-gluon Hamiltonian is build, and the structure of the glueballs 
studied is presented in Sec. IIHI The three-gluon glueball spectrum is presented with the 
two-gluon glueball spectrum from Ref. [?| and is discussed in Sec. |lVl Some concluding 
remarks are given in Sec. |V] 
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II. HAMILTON! AN 



A. Parameters 



In Ref. Q, two sets of parameters, denoted A and B, were presented for the model III 
(see Table |lj. With the set A, it is possible to obtain glueball masses in agreement with the 



results of some experimental works hJ the lowest 2 ++ state near 2 GeV, the lowest 



++ state near 1.5 GeV, and the lowest 0~ + state near 2.1 GeV. The values of a and as 
are close to the ones used in some recent baryon calculations With the set B, glueball 
masses were computed in agreement with the results of the lattice calculations of Ref. 
If the absolute glueball masses found in Ref. []| with both sets are strongly different, the 
relative spectra are nearly identical. As we use in this work a three-body generalization of 
the Hamiltonian model III of Ref. 3|, the two sets will also be considered. 

It is worth mentioning how the parameters have been determined in Ref. jjj . The mass of 
the lightest 2 ++ is nearly independent of the values of as and 7, but depends strongly on a. 
So, this last parameter has been determined with this 2 ++ state. The remaining parameters 
as and 7 have then be computed in order to reproduce the lightest ++ and 0~ + states. The 
three states 2 ++ , ++ , and _+ have been chosen because they are possible experimental 



glueball candidates 
lattice calculations 



2] 



PL lllj and because they are computed with relatively small errors in 



B. Confinement potential 

A good approximation of the confining interaction between a quark and an antiquark in a 
meson is given by the linear potential ar, where r is the distance between the two particles 
and where a is the string tension. In a baryon, lattice calculations and some theoretical 
considerations indicate that each quark generates a flux tube and that these flux tubes meet 
in a junction point -Ro which minimizes the potential energy. Following this hypothesis, the 

confinement in a baryon could be simulated by the three-body interaction 

3 

V qqq = a^2\ri-Ro\, (1) 

i=l 

For such a potential, the point Rq minimizes also the length of the three flux tubes and is 
identified with the Toricelli point [13]. 
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The energy density A c of a flux tube (string tension) can depend on the color charge c 
which generates it. Lattice calculations 14 1 and effective models of QCD predict that 
the Casimir scaling hypothesis is well verified in QCD, that is to say that the energy density 
is proportional to the value of the quadratic Casimir operator F c 2 of the colour source 

Ac = Fl a. (2) 

We have then \ q = Xq = 4<r/3 = a and \ g = 3er. In this work we will assume that the 
confinement in a three-body color singlet is given by 

3 

V ccc = aJ2 F ?\ r i-Ro\- (3) 

i=i 

This potential can be considered as the three-body generalization of the confinement used 
in Ref. hf. No constant potential is added, contrary to usual Hamiltonians in mesons and 
baryons plQ- Let us note that if the three color charges are not the same, Ro is no longer 
identified with the Toricelli point 

Interaction Q is very difficult to use in a practical calculation. A good approximation 
can be obtained for three identical color charges by replacing Rq by the center of mass 
coordinate R cm and by renormalizing the potential by a factor / which depends on the 
three-body system For three identical particles, the best value is / = 0.9515. We 

will use this approximation in the following, which seems more realistic than a confinement 
obtained by the sum of two-body forces p, t| ■ 

As already mentioned, onlyJ±Le L = states are studied in this paper. So, no spin-orbit 
correction to t 
In Refs. 



;he confinement j?J is taken into account here. 

*. BO, 3, 



the confinement potential saturates at large distances in order to sim- 
ulate the breaking of the color flux tube between gluons due to color screening effects. An 
interaction of type (J3J) seems a priori inappropriate since the potential energy can grow with- 
out limit. But the phenomenon of flux tube breaking must only contribute to the masses of 
the highest glueball states. Moreover, it has been shown that the introduction of a satura- 
tion cou 



objects 



d not be the best procedure to simulate the breaking of a string joining two colored 

3- 
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C. Dynamical constituent gluon mass 



lich can 
the 



, wind 



Within the auxiliary field formalism (also called einbein field formalism) 
be considered as an approximate way to handle semirelativistic Hamiltonians 
effective QCD Hamiltonian has a kinetic part depending on the current particle masses 
rrii and the interaction is dominated by the confinement. A state dependent constituent 
mass jj{ = (yjpf + rnf) can be defined for each particle, and all relativistic corrections (spin, 
momentum, . . . ) to the static potentials are then expanded in powers of This approach 
has been used in Ref. to build the two-gluon Hamiltonian. So, the same formalism will 
be applied also in this paper. 



Taking into account the considerations of Sec. Ill B( the simplest generalization to a three- 
gluon system of the dominant part of the model III two-gluon Hamiltonian of Ref. |7[ is 

3 3 

i=l i=l 

with the condition Y2^=i Pi = 0> since we work in the center of mass of the glueball. The 
gluons have vanishingcurrent masses and their color is such that (F?) = 3. Contrary to 
SOm e <_ wo JBflQ, o Ur Mt oaia» is a se— vistic oae la Re, & it has 
been shown that it is an important ingredient to obtain correct two-gluon glueball spectra. 

Using the technics of Ref. j^J, it is possible to obtain an analytical approximate formula 
giving the glueball mass M and the constituent gluon mass fi (the three constituent gluon 
masses are the same since the wave function is completely symmetrized, see Sec. IIII|) for the 
L = eigenstates of Hamiltonian Hq 



M ~ 6 /i with fi w 2 



7T 



+ ~, r— i /' + I 



3 x 5 x /3 9 x 5 5 / 6 V 2 



(5) 



The accuracy of these formulas, which is around 5% for the ground state (n — 0), deteriorates 
with increasing values of n. With a value of a = 4<r/3 around 0.2 GeV 2 , the smallest gluon 
constituent mass is around 600 MeV. It is then relevant to use an expansion in powers of 
1/yUo- Such a value of the gluon mass is in agreement with the values used in Refs. j^EI EJ, 
but here the constituent mass is dynamically generated. 

Instead of using the auxiliary field formalism, it is possible to consider relativistic correc- 



tions which are expanded in powers of l/E^Pi) where Ei(pi) = \/ pf + mf (see for instance 
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Ref. |2l|). But, this leads to very complicated non local potentials which are difficult to 
handle. 

D. Short-range potential 

The Hamiltonian H (jlj) gives the main features of the three-gluon glueball spectra, but 
the introduction of a short-range potential is necessary to achieve a detailed study. In 
Ref. [3], a OGE interaction between two gluons, coming from Ref. j^j], has been considered. 
It is not possible to use it directly for a three-gluon glueball because the color structure of 
the interaction is different. So, we use here the last version of a OGE interaction between 
two gluons developed specifically for three-gluon glueballs 0, Q . Its explicit form, which is 
very similar to the form of the OGE interaction for two-gluon glueballs, is given below. 

This interaction contains a tensor part and a spin-orbit part. Both are neglected in this 
paper since only L = states are studied. Moreover, it has been shown that the tensor 



interaction between two gluons is small in two-gluon glueballs j2j. 

The OGE two-gluon potential has a priori a very serious flaw: depending on the spin state, 
the short-range singular part of the potential may be attractive and lead to a Hamiltonian 
unbounded from below [6]. This problem is solved, as in Ref. |2J, by giving a finite size to 
the gluon (see Sec. IIIEJ) . 

The OGE two-gluon potential depends on the gluon constituent mass. To determine it, 
we follow the procedure proposed in Ref. Q. For a given set of quantum numbers {a}, the 
eigenstate \<p a ) of the Hamiltonian H is computed. With this state, a constituent gluon mass 
is computed fi a = (4> a \^/^\4>a)- This value of fj, a is then used in the complete Hamiltonian 




(see Sec. Ill F|) to compute its eigenstate with quantum numbers {a}. It is worth noting 
that, with this procedure, two states which differ only by the radial quantum number are 
not orthogonal since they are eigenstates of two different Hamiltonians which differ by the 
value of fi. It is shown in Ref. [9] that this problem is not serious, the overlap of these states 
being generally weak. 
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E. Gluon size 



In potential models, the gluon is considered as an effective degree of freedom with a 
constituent mass. Within this framework, it is natural to assume that a gluon is not a pure 
pointlike particle but an object dressed by a gluon and quark-antiquark pair cloud. Such an 



22j and baryon [23] sectors. As 



hypothesis for quarks leads to very good results in meson 
in Ref. we assume here a Yukawa color charge density for the gluon 

P{u) = jV , (6) 

where 7 is the gluon size parameter. The interactions between gluons are then modified by 
this density, a bare potential being transformed into a dressed one. 

The main purpose of the gluon dressing is to remove all singularities in the short-range 
part of the interaction 23]. But, for consistency, the same regularization is applied to the 
confinement potential, although no singularity is present in this case. We think that the 
definition of a gluon size, which has a clear physical meaning, is preferable to the use of a 
smearing function only for potentials with singularity 0, 0] . 

A one-body potential, like |r*j — i? cm |, is dressed by a simple convolution over the density 
of the interacting gluon and the potential 

V{r)* = J dr'V{r')p{r -r'). (7) 

A dressed two-body potential, depending on \r{ — Vj\, is obtained by a double convolution. 
This procedure is equivalent to the following calculation 

V(r)** = J dr' V{r') T(r - r') with r(«) = ^^e~ uh . (8) 
F. Total Hamiltonian 



To obtain the total Hamiltonian for three-gluon glueballs which is the simplest general- 
ization of the Hamiltonian for two-gluon glueballs from Ref. jj], we take the Hamiltonian 
-f^o given by the relation we add the OGE interactions coming from Ref. |(| (without 
spin-orbit and tensor parts); and we dress all the potentials with the gluon color density ©. 
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This gives the following Hamiltonian 

3 

H = E VP" + V OOE + ^Conf With 

'1 1 



i=l 



Voce = <* 8 ]T F, ■ F j 

i<j=i 



U{r)** 


1 

(/i 2 7 2 - l) 2 


5(r)** 


— ^ c~ r ^ 

87T7 3 


V* 

^Conf 


= 

1=1 



7T 



4 + 3^i)^r-^K-r(/ 3 +^ 



-r/-y 



+ 



r / 27(/i 2 7 2 — 1) 



with C/(r) 



-/nr 



, 1 - e" r / 7 



(9a) 

(9b) 

(9c) 
(9d) 

(9e) 



where a = 3 a/4, ^i=iP« = an d = &i + ^i- /5 = +1 (—1) for a gluon pair in color 
octet antisymmetrical (symmetrical) state. The constituent state-dependent gluon mass \i 
is computed in advance with a solution of the Hamiltonian Hq. 



III. WAVE FUNCTIONS 

A gluon is a I(J P ) = 0(1 _ ) color octet state. Two different three-gluon color singlet states 
exist p], which are completely symmetrical or completely antisymmetrical ((Fi ■ Fj) = 3 
for such states). The total isospin state of a glueball is an isosinglet and is completely 
symmetrical. Different total spin states are allowed with different properties of symmetry. 
They are presented in Table |H] As gluons are bosons, the total wave function must be 
completely symmetrical. Its parity is the opposite of the spatial parity, and its C-parity 
is positive for color antisymmetrical state and negative for color symmetrical state. Let us 
note that a two-gluon glueball has always a positive C-parity. 

In this work, we will mainly consider glueballs with the lowest masses. These states are 
characterized by a vanishing total orbital angular momentum L = and by a spatial wave 
function completely symmetrical with a positive parity. This immediately implies that the 
lowest glueballs are states with J PC equal to CT + , 1 , and 3 In order to reach 

a good accuracy, the trial spatial wave functions are expanded in large gaussian function 
bases 2||. With more than 10 gaussian functions for each color/isospin/spin channels, we 
have checked that the numerical errors on masses presented are around or less than 1 MeV. 

Using the value of a from models A and B, the L = ground state masses of the 
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Hamiltonian H (J3J) are presented in Table IIHI for all possible J PC quantum numbers. The 
0~ + , 1 , and 3 glueballs have clearly the lowest masses. In this table, they are degenerate 
since the Hamiltonian Hq is spin- independent. For each state, the corresponding constituent 
gluon masses /j,q is indicated. It is used to define the complete Hamiltonian H (0). 



IV. RESULTS 



We present here the three-gluon glueball masses obtained with the complete Hamiltonian 
H (jOJ) together with the two- gluon glueball masses computed in Ref. Q| (see Table ITVL These 
masses are compared with the results obtained by the lattice calculations of Ref. 2l. This 

n q 

work is an update of a previous study So, a state not computed in Ref. |2| but presented 
in Ref. is also considered here. As it can be seen on Fig. Q with the set B of parameters, 
our masses are in quite good agreement with the results of the lattice calculations, except 
for one exception discussed below. Unfortunately, the masses predicted for the ++ , 2 ++ , 
and H two-gluon glueballs are larger than some possible experimental candidates |lfl . 
A best agreement with these data can be achieved with the set A of parameters But the 



situation is not simple, since recent works 1261 1271 suggest that the new observed resonance 
/o(1810) reported by the BES collaboration [28] could be a ++ glueball. Such a data is 
then compatible with the set B of parameters. 

As the lattice results can suffer from quite large error of normalization, it is then inter- 
esting to present a relative spectra. In Fig. the ratios of two- and three-gluon glueball 
masses on the lowest 2 ++ mass are presented. At this scale, results coming from sets A and 
B of parameters cannot be distinguished. 

Let us first point that the 1 and 3 states occur only in three-gluon glueballs, whereas 
the 0~ + state is available for two- and three-gluon systems. The mixing between these two 
channels is ignored here. Our results share some similarities with other potential models. 
For instance, in Ref. j^, the three lowest three-gluon glueballs are those with J PC equal to 
h , 1 , and 3 , but they are found around 2400 MeV and the mass splitting between 
theses states is around 50 MeV. In Ref. Q|, the same three lowest states are found: the 1 
glueball is predicted in the range 3500-3700 MeV and the mass difference between these 
three sates is around 100 MeV. 

A detailed comparison between the results of our potential model and the lattice calcula- 
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tions of Ref. |l| have been performed in Ref. [t| for the two-gluon glueballs. The conclusions 
of this work are not changed by the new lattice predictions obtained in Ref. jsj. Let us 
then focus our attention on the three-gluon glueball spectrum obtained with the set B of 
parameters. Without any new parameters, the 1 and 3 glueballs are in good agree- 
ment with the lattice predictions of Ref. |2j (see Table ITVl and Fig. QJ). The agreement is 
slightly less good for the relative spectra (see Fig. |2J) because the error is smaller for lattice 
mass ratios. Our results also suggest the existence of a three-gluon 0~ + glueball near these 
two last states. The others 0~ + states already computed by the lattice calculations can be 
identified as two-gluon systems by our model. 

The lattice results predict a 2 state at 4010 MeV near the 1 and 3 states. With 
our Hamiltonian, a mass more than 1 GeV above is computed. It is unavoidable in our 
model, since a spin 2 function has a mixed symmetry which implies a mixed symmetry for 
the space function and then a greater mass for the corresponding glueball, in agreement 
with the results of Refs. 0, |fj]. It has been checked that the spatial wave function of the 
2 state is dominated by a configuration in which each internal variable is characterized 
by one unit of angular momentum. 

We have no firm explanation for such a discrepancy. This problem could arise because the 
gluon has a constituent mass within our formalism. So, it possesses a spin as any massive 
particle, that is to say three states of polarization. In lattice calculations, the gluon is a 
massless particle with a definite helicity and then only two states of polarization. The same 
phenomenon could be at the origin of another puzzling-to some extent opposite-situation 
in the two-gluon glueball sector: the existence of 1 PC states with low masses in our model 
[3], which are actually not predicted below 4 GeV by lattice calculations [lj. Further studies 
are necessary to clarify the situation. 

The lattice results predict also several three-gluon J + ~ states in the range 2980-4780 MeV. 
These states with positive parity have a negative spatial parity and then are expected to 
have masses larger than the lowest three-gluon J states. This is manifestly not the case 
in these lattice calculations. In Ref. [3], it is shown that the structure of some two-gluon 
glueballs can be explained with our potential model by the existence of strong spin-orbit 
forces coming from the OGE interaction and the confinement. It is possible that similar 
forces act in three-gluon system to lower the masses of some of these states. 



10 



V. CONCLUSION 



The masses of pure three-gluon glueballs have been studied with the generalization of 
a semhelativistic potential m0 del Q which give, pure two-glnon gluebah spectra in good 
agreement with lattice calculations [JJ, |2j|. The short-range part of the potential is the 
sum of two-body OGE interactions. For the confinement, a potential simulating a genuine 
Y-j unction is used and it is assumed that the Casimir scaling hypothesis is well verified. 
The gluon is massless but the OGE interaction is expressed in terms of a state dependent 
constituent mass. The Hamiltonian depends only on 3 parameters fixed in Ref. [2J: the strong 
coupling constant, the string tension, and a gluon size. All masses have been accurately 
computed with an expansion of trial states in gaussian functions |25j ]. 

In this work, only the negative (natural) parity L = three-gluon glueballs are studied. 
The masses of the lowest 1 and 3 glueballs predicted by our potential model are in 
agreement with the results of a recent lattice calculations Q], but the lowest 2 state is 
found at higher energy in ag r ee m ent with other potent n^s g Q . The origin of such a 
discrepancy between both approaches is not known. It could be due to the fact that gluons 
have constituent non vanishing masses in our approach. They are then characterized by a 
spin, and not by a helicity as it could be expected for particles with a vanishing current 
mass. 

Other positive parity three-gluon glueballs predicted by lattice calculations are not con- 
sidered here j^J. We think that their properties could be explained by the action of strong 
spin-orbit forces, similar to those present into two-gluon glueballs. To take into account these 
interaction, it is necessary to consider the spin-orbit forces arising from OGE interactions 
but also from the Y-junction confinement. Such a work is in progress. 
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TABLE I: Parameters for models A and B (<r = 3 a/4). For both models, the ; 


^luon current mass 


is zero and / = 0.9515. 






Model A 




Model B 


a 0.16 GeV 2 




0.21 GeV 2 


a s 0.40 




0.50 


7 0.504 GeV" 1 




0.495 GeV" 1 


TABLE II: Characteristics of three-gluon spin functions with total spin S, intermediate couplings 


5i n t, and symmetry properties which can be obtained by coupling 


(A: Antisymmetrical, S: Sym- 


metrical, MS: Mixed symmetry). 






S S[ n t 




Symmetry 







1 A 


1 0, 1, 2 




1 S, 2 MS 


2 1, 2 




2 MS 


3 2 




1 S 


TABLE III: L = ground state masses Mq of the Hamiltonian H 


(HJ) as a function of the J PC 


quantum numbers. The corresponding constituent gluon masses fiQ 


are also given. Values in MeV 


are computed with the value of a from models A/B. The lowest masses are printed in italic. 


J PC M /xo J PC 


M 


Mo 


0— 5574/6385 929/1064 _+ 


3211 /3679 


535/613 


1— 3211/3679 535/613 1~+ 


4156/4761 


693/794 


2— 4156/4761 693/794 2~+ 


4156/4761 


693/794 


3^ 3211/3679 535/613 3~+ 


5574/6385 


929/1064 
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TABLE IV: Glueball masses in MeV and (glueball mass ratios normalized to lightest 2 ++ ). The 
two-gluon masses are taken from Ref. j?|. The error bars for lattice mass ratios are computed 
without the normalization error on the masses. The lightest ++ , 2 ++ , and ^ states are taken 
as inputs to fix the parameters. The first column indicates the valence gluon content as predicted 
by our model. 





Lattice [Ref.] 






Model A 


Model B 


gg o++ 


1710±50±80 


(0.72 ±0.03) 


\2) 


1604 


(0.78) 


1855 (0.78) 




2670±180±130 


(1.12 ±0.09) 


n] 


2592 


(1.26) 


2992 (1.26) 


2 ++ 


2390±30±120 


(1.00 ± 0.03) 


\2] 


2051 


(1.00) 


2384 (1.00) 


o-+ 


2560±35±120 


(1.07 ±0.03) 


12} 


2172 


(1.06) 


2492 (1.05) 




3640±60±180 


(1.52 ±0.04) 


n] 


3228 


(1.57) 


3714 (1.56) 


2-+ 


3040±40±150 


(1.27 ±0.03) 


12} 


ZD 1 O 


(~\ 9^ 


2984 (1.25) 




3890±40±190 


(1.63 ±0.04) 


n] 


3345 


(1.63) 


3862 (1.62) 


3++ 


3670±50±180 


(1.54 ±0.04) 


\2} 


3132 


(1.53) 


3611 (1.51) 


ggg 1~~ 


3830±40±190 


(1.60 ±0.04) 


12} 


3433 


(1.67) 


3999 (1.68) 


2™ 


4010±45±200 


(1.68 ±0.04) 


12} 


4422 


(2.16) 


5133 (2.15) 


3— 


4200±45±200 


(1.76 ±0.04) 


\2} 


3569 


(1.74) 


4167 (1.75) 


o-+ 








3688 


(1.80) 


4325 (1.81) 
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O Lattice (2006) 

• Lattice (2006) - reference 

□ Lattice (1999) 

<S> Model B 



2 2 3 



1~ 2 3 



FIG. 1: Glueball masses given in MeV. Dotted diamonds: Results from model B_(two-gluon masses 
are taken from Ref. 



|); Black and white circles: Lattice results from Ref. 







White squares: 



Lattice results from Ref. |2(]. Black circles indicate the reference states taken as inputs to fix the 
parameters. The error bars for lattice results are computed by summing the two uncertainties (see 
Table E). 
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O Lattice (2006) 

# Lattice (2006) - reference 

□ Lattice (1999) 

<3> Model B 



5» 



♦ 5 



2 2 3 



! 2" 3" (T + 



FIG. 2: Glueball mass ratios normalized to the lightest 2 ++ state (see Table HVj) . Dotted diamonds: 
Results from model B (two-gluon masses are taken from Ref. 7]); Black and white circles: Lattice 



D 



results from Ref. lj; White squares: Lattice results from Ref. |2|. Black circles indicate the reference 
states taken as inputs to fix the parameters. The error bars for lattice results are computed without 
the normalization error on glueball masses. Results from models A and B cannot be distinguished 
on this graphic. 
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